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Low-capacitance Josephson junction systems as well as coupled quantum dots, in a parameter 
range where single charges can be controlled, provide physical realizations of quantum bits, discussed 
in connection with quantum computing. The necessary manipulation of the quantum states can be 
controlled by applied gate voltages. In addition, the state of the system has to be read out. Here 
we suggest to measure the quantum state by coupling a single-electron transistor to the q-bit. As 
long as no transport voltage is applied, the transistor influences the quantum dynamics of the q-bit 
only weakly. We have analyzed the time evolution of the density matrix of the transistor and q-bit 
when a voltage is turned on. For values of the capacitances and temperatures which can be realized 
£N) ' by modern nano-techniques the process constitutes a quantum measurement process. 
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I. INTRODUCTION 

Recent proposals&B suggested to use nanoscale devices, such as low-capacitance Josephson junctions or coupled 
quantum dots as quantum bits (q-bitsi, which are the basic elements of quantum computers. The two logical states are 
different charge states of the systemlini Applied gate voltages allow the necessary controlled manipulations (single-bit 
and two-bit operations) of the quantum states. In addition to these manipulations, a read-out device is required to 
perform quantum measurements of the resulting state of the q-bit. We suggest to use single-electron transistors for 
£h ' this purpose. 

The requirements to perform, on one hand, quantum manipulations and, on the other hand, a quantum measurement 
appear to contradict each other. During the manipulations the dephasing should be minimized, while a quantum 
measurement should dephase the state of the q-bit as fast as possible. The option to couple the measuring device to 
the q-bit only when needed is hard to achieve in mesoscopic systems. The alternative, which we discuss here, is to 
keep the measuring device permanently coupled to the q-bit in a state of equilibrium during the quantum operations. 
The measurement is performed by driving the measuring device out of equilibrium, in a way which dephases the 
quantum £tate of the q-bit. Similar nonequilibrium dephasing processes have recently been considered by a number 
jy-^ , of authoraJH. 

' For definiteness we discuss in this paper the measurement process performed by a single-electron tunneling (SET) 
t-H \ transistor coupled capacitively to a Josephson junction q-bit; however, this type of measurements may be performed 
t-H ' for any quantum system with two different charge states. We describe the measuring process by considering the time- 
evolution of the density matrix of the coupled system. We show that the process is characterized by three different 
time scales: the dephasing time, the time of measurement, which may be longer than the dephasing time, and the 
mixing time, i. e. the time after which all the information about the initial quantum state is lost due to the transitions 
induced by the measurement. Thus, we arrive at a new criterion for a "good" quantum measurement: the mixing 
time should be longer than the time of measurement. 



' 

£ 1 II. THE QUANTUM SYSTEM AND THE MEASURING DEVICE 

o : 

The system is shown in Fig. [lj The two superconducting islands in the upper part are the realization of a q-bit. Its 
state is characterized by a discrete variable, n, the number of extra Cooper pairs on the lower superconducting island. 
The lower part (a normal island between two normal leads) stands for a SET transistor, which is coupled capacitively 
to the q-bit. It's charging state is characterized by the extra charge on the middle island, eN. A similar setup has 
recently been studied in the experiments of Refs. [ p|jic|l with the purpose to demonstrate that the ground state of a 
single Cooper pair box is a coherent superposition of different charge states. We discuss the relation of our proposal 
to these experiments below. 

As shown earlieiu the quantum operations with the q-bit are performed by controlling the applied gate voltage V q b- 
At this stage the transport voltage V across the SET transistor is kept zero. Therefore no dissipative currents flow 
in the system, and dephasing effects due to the transistor are minimized. To perform a measurement one applies a 
transport voltage V. The resulting normal current through the transistor depends on the charge configuration of the 
q-bit, since different charge states induce different voltages on the middle island of the SET transistor. In order to 
check whether the dissipative current through the SET transistor contains information about the quantum state of 
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the q-bit, we have to discuss various noise factors (shot noise) and the measurement induced transitions between the 
states of the q-bit. It turns out that for suitable parameters, which can be realized experimentally, the dephasing by 
the passive SET transistor is weak. When the transport voltage is turned on the dephasing is fast, and the current 
through the transistor - after a transient period - provides a measure of the state of the q-bit. At still longer times 
the complicated dynamics of the composite system destroys the information of the quantum state to be measured. 

The Hamiltonian of the composite system consists of three main parts: the charging energy, the terms describing the 
microscopic degrees of freedom of the metal islands and electrodes, and the tunneling terms, including the Josephson 
coupling. The charging term is a quadratic form in the variables n and N: 

^charge = E qh n 2 + E sct N 2 + E int nN + 

2enV n + eNV N + const. . (1) 

The charging energy scales E qbl E sct and i?i n t are determined by the capacitances between all the islands. Similarly, 
the effective gate voltages V n and Vn depend in general on all three voltages Vqb, V g and V, but for a symmetric bias 
(see Fig. ||), V n and Vn are controlled only by the two gate voltages, Vqb and V g . 

The microscopic terms , and H\ describe noninteracting electrons in the two leads and on the middle island 
of the SET transistor, respectively: 

Hr = J2^ c kA* (r = L,R,I), (2) 

ha 

The index a labels the transverse channels including the spin, while k labels the wave vector within one channel. 
Similar terms exist for the two islands of the q-bit. Here we use the "macroscopic" description of the superconductors, 
assuming that the microscopic degrees of freedom have already been integrated outO. 

The tunneling terms include the Josephson coupling Hj = —Ej cosO, which describes the transfer of Cooper pairs 
between the two islands of the q-bit (e ie \n) = \n + 1)), and the normal tunneling Hamiltonian for the SET transistor: 

Ht = Y, TwAUl'o^ + E ^V.4U v e-^ + h.c. . (3) 

kk' a k' k" a 

Apart from the microscopic degrees of freedom, (||) contains two "macroscopic" operators: e ±l ^ and e ±l ^ . The first 
one describes changes of the charge on the transistor island due to the tunneling: e^lN) — \N + 1). It may be treated 
as an independent degree of freedom if the total number of electrons on the island is large. We include one more 
operator e ±l ^ which describes the changes of the charge in the right lead. It acts on m, the number of electrons 
which have tunneled through the SET transistor, e"^|m) = |m + 1). Since the chemical potential of the right lead 
is controlled, m does not appear in any charging part of the Hamiltonian. However, e ±4 ^ allows us to keep track of 
the number of electrons which have passed through the SET transistor, which is related to the current through the 
device. 

We define the q-bit 's Hamiltonian as the part of the total one which governs the q-bit 's dynamics in equilibrium 
(N = 0): 

H qh = E qh {n - Q qh f - Ej cos 9 . (4) 

Here Q q b = — 2eV n /2E q t, is the q-bit's gate charge, measured in units of 2e. We concentrate on the values of Q q b m 
an interval around the degeneracy point Q q b = 1/2, so that only the low energy charge states n = and n = 1 are 
relevant. These states, however, are not appropriate logical states of the q-bit since they are not the eigenstates of the 
Hamiltonian (||). We diagonalize (||) in the two charge states subspace for a fixed value of Q q b (which is kept constant 
between the quantum manipulations and during the measurement) and denote the corresponding logical states |0) 
and |1). In the new basis, up to a constant, H q h = — (l/2)AEa z , where a z is the Pauli matrix and 



AE= ^/[E qh {l-2Q qh )} 2 + E] . (5) 

The price which we pay for this simplification is that the number operator n, which appears in the mixed term of (Q), 
becomes non-diagonal: 

n= - — -cos rj a z - -sm r\a x , (6) 

with mixing angle rj given by tan 77 = Ej/E q b(l — 2Q q b). In the quantum regime, which we are considering here, 
E q b ^> Ej and, therefore, one can choose Q q b so that tan r] -C 1. 
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The interaction Hamiltonian (part of the mixed term in (Q)) now becomes 

H int = - -E int N (cos r] (j z + sin na x ) , (7) 

while the rest of the mixed term (E- m tN/2) as well as all other remaining terms are collected in the Hamiltonian of 
the SET transistor: 

H sct = E sct (N - Q sct ) 2 + H L + H R + H I + H T . (8) 

The transistor's gate charge (measured in the units of e) became Q so t = — (eVN+Ei nt /2)/2E sct . The total Hamiltonian 
reads H = H q \, + H sct + H lnt . One should understand, however, that the division chosen is rather arbitrary. The 
terms -ff q b and H sct would not describe the q-bit and the SET transistor if they were decoupled. 

III. QUANTITATIVE DESCRIPTION OF THE MEASUREMENT, TIME EVOLUTION OF THE 

REDUCED DENSITY MATRIX 

The total system is described by a reduced density matrix o(t) = TrL,R,i{/9(t)}, where the trace is taken over the 
microscopic states of the leads and of the island. In general, the density matrix a(i,j;N,N';m,m') is a matrix m. 

which stand, for the quantum states of the q-bit (|0) or |1)), in N, and in m. However, as has been shown intHI 
a closed set of equations describing the time evolution of the system can be derived where the off-diagonal elements 
in N have been eliminated. The same is true for the off-diagonal elements in m. Therefore, we need to consider only 
the following elements of the density matrix crf- ra = cr(i, j; N, N;m,m). We assume now that at time t — 0, when 
the q-bit is prepared in the quantum state a|0) + 6|1) as a result of previous quantum manipulations, we switch on a 
transport voltage to the SET transistor. To proceed we can further reduce the density matrix in two different ways 
to obtain dual descriptions of the measuring process. 

The first widely used procedural] is to trace over N and m. This yields a reduced density matrix of the q-bit 
a i,j = m a ij m ■ Assuming that at t = it is in the state 

.ww- (Kg:). < 9 > 

the questions arc how fast the off-diagonal elements of a^j vanish (dephasing), and how fast the diagonal elements 
change their original values (for instance due to transitions induced by the measurement). This description is enough 
when one is interested in the quantum pxaDerties of the measured system only (q-bit in our case) and the measuring 
device is used as a source of dephasingd'EHj. It does not tell us much, however, about the quantity measured in an 
experiment, namely the current flowing trough the SET transistor. 

The second procedure is to evaluate the probability distribution of the number of electrons m which have tunneled 
trough the SET transistor during time t: 

P(m,i) = $>£ m (t). (10) 

N,i 

This quantity gives a complete description of the measurement. At t = no electrons have tunneled, so P(m, 0) = 5 m $. 
Then this delta-peak starts to shift in positive m direction and, at the same time, it widens due to shot noise. Since 
two states of the q-bit correspond to different conductivities (and shift velocities in m space) , one may hope that after 
some time the peak splits into two. If after sufficient separation of the two peaks their weights (integrals) are still 
close to | <z | 2 and |6| 2 , a good quantum measurement has been performed. Unfortunately, there exist further processes 
which destroy this idealized picture. After a long time the two peaks transform into a broad plateau, since transitions 
between the q-bit's states are induced by the measurement. Therefore, one should find an optimum time for the 
measurement, so that, on one hand, the two peaks are separate and, on the other hand, the induced transitions have 
not yet happened. In order to describe this we have to analyze the time evolution of the reduced density matrix 
quantitatively. 

IV. DERIVATION OF THE MASTER EQUATION 

The Blftch-type or master equations with coherent terms have only recently been analyzed in the condensed matter 
physical3i3. In Rcf. [ [l2| a diagrammatic technique has been developed which provides a formally exact master 
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equation as an expansion in the tunneling strength. Only the tunneling term Ht is considered a perturbation, while 
all other terms constitute the zeroth order Hamiltonian Hq = H — H^, which is treated exactly. The master equation 
reads 



da(t) 
dt 



~[a(t),H ] = J dt'X(t-t')a(t') , 



(11) 



where the_matrix elements of T,{t',t) can be calculated diagrammatically using the real-time Keldysh contour 
techniquet3. The simplest diagram describing the tunneling through the left junction in first order perturbation 
theory (secuiential tunneling) is shown in Fig. @. The dashed lines crossing the diagram contribute the following factor 
to the rateB: 



(JE_'* 2 p±V'l(*-*') 



sinh 



Mt-t'±iS) 



(12) 



where oil = 7i/(47r 2 e 2 i?T L ), Ml is the electro-chemical potential of the left lead, and 5 = l/u) c is the inverse frequency 
cut-off. The sign of the id term depends on the time-direction of the dashed lincO. It is minus if the direction of the 
line with respect to the Keldysh contour coincides with its direction with respect to the absolute time (from left to 
right), and plus otherwise. For example the right part of Fig. || should carry a minus sign, while the left part carries 
a plus sign. For the sign in front of ifx^t — t') the rule is as follows: minus, if the line goes forward with respect to 
the absolute time, and plus otherwise. i_. 

For a single SET transistor the horizontal lines correspond to trivial exponential factorsEj e lEt . In our case, however, 
we have to account for the nontrivial time evolution of the q-bit. Therefore the upper line in the left part of Fig. ^ 
corresponds to (N — 1, j\e~ %H °( t ~ t ^\N — while the lower line corresponds to (N, * | e*^ (* — * )|JV,i'). To calculate 

these matrix elements we diagonalize iJ q b + -ffint for each value of N. The eigenenergies are 



4? = T^A^ + ^intiVcosry) 2 + (£ int 7Vsinr/) 2 



(13) 



and the mixing angles en (analogous to 77) are given by tane^ = £ , i nt -/Vsin?7/(A_E-|-£ , i n t J ?V cos 77). The matrix elements 
(propagators) read: 



(N,0\e 


-* H ° At \N, 0) 


(N,l\e 




(N,l\e 





,2 CN -iE»M 



iE" At 



sm 



2 £W - iE ? Au -iE^At 



-iE. v AtN 



sine N (e ~ - < ~ ) < 



Ei At \ -iE'Z'At 



(14) 



where = E sct (N - Q sct ) 2 . 

We analyze now the rates in Fig. || for different choices of q-bit's indices in the regime AE 3> E- mt , Ej. There the 
mixing angles are small, ejy oc N E- lnt E j / (AE) 2 , for all relevant values of N. Hence, we keep only terms linear in e^r. 
The simplest transition (i 1 = 0, j' = 0, N — 1, m) — * (i = 0, j = 0, N, m) is described by 



l,m,0; A/", m,0 / a ,\ 
l,m,0;iV,m,0^ t J 



-a L (^) 2 e 



iEAt 



sinh 



hr) 



(At + iS) 



c.c. 



(15) 



where E stands here for ML + (E^~ 1] - Ej£)) + (E^ 1 - Eg). 

The form of the master equation ( |ll"| ) suggests the use of a Laplace transformation, after which the last term in 
( |j"l| ) becomes £(s)er(s). We Laplace transform ( |l5| ) in the regime s <C E, i.e. we assume the density matrix a to 
change slowly on the time scale given by h/E. This assumption should be verified later for self-consistency. At zero 
temperature (j3 — > 00) and for S — > we obtain: 



-,(1) iV— l,m,0;N,m,0 



N-l, 



.,0:JV,t 



(s) = 2a L Re 



(s + iE)e* 5 ^ a+iMr >E 1 [iS(a + iE)} 



2Tra L EO(E) - 2a L s(l + 7 + M\SE\)) 



(16) 



where J5i[...] is the exponential integral and 7 ~ 0.58 is Euler's constant. Denoting the diverging factor [l+7+ln(|5£J|)] 
by D(E) and performing the inverse Laplace transform we arrive at: 
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n-Mn'X'°o(^) « 2 7 ra LJ Be(^) ( 5(Ai + 0) - 2a L D(E)5'(At + 0) . (17) 

(Note that (O) is equivalent to ( |lE|) only as a kernel in the convolution ([ll]) when applied to slowly changing matrix 
elements of a). The first term of ( |17| ) is the usual Golden Rule tunneling rate corrected with respect to the additional 
charging energy corresponding to the quantum state |0) of the q-bit, E^ -1 — Eq . The second (diverging) part of ( p"7| ) 
produces a term proportional to ^ a ~ 1,m . One can take this term to the LHS of ( pd| ) so that the time derivative in 
the LHS will look like Tt [a Q 

o 1,m ]- We analyze all possible choices of the q-bit's indices in Fig. || and 
arrive at the conclusion that the diverging terms have always the same structure as the coherent terms in the LHS 



of (11). Moreover, if we neglect some energy corrections of order i?i n t, we may incorporate all of these terms to the 



LHS of (|1 1|) , so that the master equation reads: 



(1 + a L A + a R B) 



= T*(t) , (18) 



where A and B are three-diagonal matrices in the N and m spaces, composed of the diverging factors of the type of 
D(E), while T is the regular local part of £(i — t'). 

We expect that without the approximation of energies in the diverging terms the structure of ( |l8| ) would be the 
same, with A and B being more complicated matrices, which would include some mixing in the space of the q-bit's 
states. Finally, we note that for any physically reasonable choice of the cut-off S, the logarithmically divergent factors 
in the matrices A and B are of order one, and, therefore, the mixing corrections to the unit matrix in the LHS of (|lg| ) 
are small. We multiply the master equation ( p"8| ) by (1 + a^A + a^B)^ 1 « (1 — a^A — orB) from the left, so that the 
mixing corrections move to the RHS. Since T is linear in «l and «r, the mixing corrections are quadratic. We drop 
them in the framework of the first order perturbation theory. The master equation to be analyzed thus becomes: 

^--L[a(t),H ]=Ta(t). (19) 

If the applied voltage is not too high (the exact criterion to be specified) we may consider only two charge states 
of the SET transistor, N = 0,1. We perform a Fourier transform in m space cr^-(fc) = J2 m a i j m e lkm . To shorten 
formulas we introduce A N = a^ (k) , B N = erf x (k) , C N = £ Re <j*[ m e tkm , and D N = £ Im a^f 1 e lkm . This enables 

m m 

us to rewrite (|19|) as: 

A = -T ho A° + Tn^A 1 - lo l C° - u K e %k C x (20) 

A 1 = T Lo A° - r Ro A 1 - SID 1 + lo l C° + conC 1 (21) 

B° = -r Ll B a + L Rl e ik B l - w L C° - lu^C 1 (22) 

b 1 = r Ll b" - r Rl b 1 + nD 1 +uj- L c° + ^rC 1 (23) 



C° = -AE°D° - T L C° + Tne^C 1 - ^ (A - B°) - ^"(A 1 - B 1 ) (24) 
C 1 = -AE l D l + T L C° - TnC 1 + ^ (A - B°) + ^-(A 1 - B 1 ) (25) 



— -Uli 1/ TIL 1 - 

D° = AE°C° -T L D° + T R e lk D 1 (26) 
D 1 = AE 1 ^ + T L D° - TrD 1 + ^(A 1 - B 1 ) . (27) 

Here AS 0,1 = E^' 1 — E®' 1 are the energy differences between the q-bit's states for N = and N — 1 respectively, 
and Q = E- ln t sin 77 is the coefficient in the mixing term in H lnt for N — 1 (see ((?])). The terms proportional to AE ' 1 
and fl originate from the coherent part of (p^). The tunneling rates which appear in the four last equations for the 
off-diagonal elements ( p4| - ^7| ) are given by 

T L = 27ra L [/i L - E sct (l - 2Q set )} 

T R = 27ra R [-/i R + E sct {l - 2Q set )] (28) 



The rates 



r Lo>1 = t l ± Ar L 

r Ro , = T R ± AT R (29) 
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appearing in the equations for the diagonal elements (|20||23|) are corrected due to the charging energy induced by the 
q-bit: 

Ar L = 2 TOL (^-^ 1 ) 

Ar R = -2 7 ra R (^-^ 1 ) . (30) 

These correction are, actually, responsible for the separation of the peaks. In the regime tan 77 -C 1, which we assume 
here, |ATl,r| ~ 27rctL,R-E'int- The rest are small mixing terms, o->l,r = i"Q!l,r sineiAi? 1 , which also originate from 
diagrams of the type in Fig. ||. Note, that we assume that only the two rates given in ( p8| ) are nonzero (two-state 
approximation). Moreover we assume that the q-bit's charging energy corrections can at most change these two rates, 
but they can not switch on any other rate or switch off one of the two in (E8h . 



V. QUALITATIVE ANALYSIS OF THE MASTER EQUATION 



First, we analyze the system (|20j|27]) qualitatively. Imagine that we can "switch off" the Josephson coupling during 
the measurement. Then all the mixing terms in (|20|-p7|), i.e. those proportional to and ^l.r disappear, and the 
system factorizes into three independent groups. The first one (p0||2l|), 



A 



-T Lo A- ti Rl 
A 1 =T Ul A° -T^A 1 , 

has plain wave solutions (oc e luJt ). The standard analysis gives for eigenvalues: 



(31) 



W12 



r(r Lo +r Ro Wi± 



4r ( e 



ik 



Ro 



(32) 



where 



r n = 



_ r Lo r Ro 



r Lo + r 



Ro 



(33) 



is the total transport rate corresponding to the q-bit is in the state |0). When k is small u>\ ~ iij^h + Tro); while 
u>2 ~ Tofc + iTo(l + „ 21 ^p — )k 2 . Since wi is a large imaginary number, already after a short time, l/lwil, only the 

second eigenvector (A 1 /A° = Fl /F Ro ) survives. This eigenvector multiplies a wave packet propagating with the 
group velocity r . The wave packet widens due to shot noise of the single electron tunneling, it's width being given 
by v^Fot (the second imaginary term in the expression for 0J2). 



Tl 1 Tr 1 



and 



Analogously the second group of equations (|22| , |23D gives a wave packet with the group velocity Ti = (p L +Fn 

the width w y/Tit. The two peaks correspond to the q-bit in the states |0) and |1), respectively. They separate when 
their width is smaller then the distance between their centers y/Toi + \fTii < \Tq — Ti\t. After this time, 

(34) 



tms — IvTo — v ri 



which we denote as the measurement time, the process can constitute a quantum measurement. Similar expressions 
have been obtained in Refs. [ H-f?J, where they have been denoted as dephasing time. 

To get a clue for the dephasing we analyze the third group of equations ( p4[]27| ) at k = (the trace over m is 
equivalent to k — 0). These equations may be combined into two complex ones: 



(35) 



I AE 1 — AE a \ m Ei nt <C (Tl + Tr) the imaginary parts of the eigenvalues are 
In the opposite limit dE (Tl + Tr) the imaginary parts are Imwi ~ Tl 



The standard analysis shows that if dE 
Imwi w (Tl + r fl ) and Imw 2 ~ 4 ( r ^+r fl ) • 
and Imw2 ~ Tr. The first limit is physically more relevant (we have assumed parameters in this regime), although 
the second one is also possible if the tunneling is too weak or the coupling between the q-bit ant the SET transistor 
is too strong. In both limits the dephasing time, which is defined as the the longer of the two times, 



G 



= max{ [Im u>x] \ [Imu^] *} (36) 
is parametrically different from the measurement time j34|) . In the first limit, dE <C (r^ + Fr), it is 

4(r L + v R ) 

T 4> = ^2 K a L,R (37) 

while < ms cx «l p. One can check that in the whole range of validity of our approach the measurement time exceeds 
the dcphasing time, i ms > r^. This is consistent with the fact that a "good" quantum measurement should completely 
dephase a quantum state. In Refs. [^-0, where different systems have been discussed, the expressions for the resulting 
dephasing time were given by expressions similar to ([mJ), thus t$ — t ms . 

In our example the dephasing time is shorter than the measurement time. The reason for this is, probably, the 
presence of the additional uncontrolled environment provided by the middle island of the SET transistor. The 
transport of electrons occurs via a real state of the island N = 1 . In different transitions the island may be left in 
different microscopic states, even though the same number of electrons have passed. To put it in the language of 
Ref. [ 0, the initial state of the system (a\Q) +b\l)) \x) \m — 0) evolves into a|0) \xo) \ m o) + b\l) \xi) \mi)> where \x) 
stands for the quantum state of the uncontrolled environment. One may imagine a situation when mo = mi, but |xo) 
and \xi) are orthogonal. In this situation the dephasing has occured but no measurement has been performed. 

The additional environment plays, actually, a positive role, i. e. it helps us to perform a quantum measurement, 
provided it dephases the state of the q-bit only when the system is driven out of equilibrium. This is because the 
dephasing suppresses the transitions between the states of the q-bit (Zeno effect). 



VI. THE MIXING TIME 

Finally, we analyze what happens if we take into account the mixing terms in the system ( ^p|^7|). We assume 
k = and investigate the eigenvalues of the eight by eight matrix formed by the coefficients of (|20|-[27|). Note that 
in the discussion above we have calculated all the eight eigenvalues for _Ej = (the two eigenvalues of the complex 
system ( |35| ) are doubled when one considers it as a system of four real equations) . In the diagonal part there were 
two zeros, which corresponded to two conserved quantities, A + A 1 = ero,o an d B a + B 1 = o\.\. Six other eigenvalues 
were large compared to the amplitudes of the mixing terms. It is clear, that switching on the mixing changes only 
slightly the values of the six large eigenvalues. Moreover, one of the eigenvalues is always zero. This corresponds to 
the conservation of the total trace A a + A 1 + B a + B 1 = 1. The last (8th) eigenvalue acquires now a small imaginary 
part and this gives the time scale of the mixing between the two states of the q-bit. 

We do not have an analytical expression for the mixing time, but we can estimate it for a concrete physical 



situation. At the degeneracy point, we have IY = Tr, and the corrections to the rates (30) cancel each other, thus, 
no measurement is performed. Therefore, we choose Q se t far enough from the degeneracy point, which is Q sc t = 1/2, 
so that Fl < Tr and the Coulomb blockade energy Eqb = E sct (l — 2Q sct ) is of the order of E set . To satisfy the 
conditions for the Golden Rule (see (|l|) and the discussion thereafter) we assume the chemical potential of the left 
lead /iL = V/2 to exceed the Coulomb blockade energy by an amount of the order of Eqb °c E sct and assume Eqb to 
be the largest energy scale of the system: Ecb A£\ The transport voltage should not, however, exceed the limit, 
after which the third charge state of the SET transistor N = — 1 becomes involved. Thus V/2 < E sct (l + 2Q sct ) and 
Q S et should be chosen far enough from zero as well. In this regime we estimate the mixing time as 

(AE) 4 



where «l = c*r = ot. The measurement time in the same regime is given approximately by 



t~l cx 2na§± . (39) 

"set 



The exact values of Q se t and V would determine the numerical coefficients in front of ( p8[ ) and (|39|). Thus, £ m s/imix oc 
Ej E^ et I (AE) 4 . One recognizes two competing ratios here: Ej/AE, which is small, and E sct /AE, which is large. The 
condition i ms /i m ix ^ lj thus, imposes an additional restriction on the parameters of the system. 
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VII. DISCUSSION 



To show that all the conditions assumed in this paper are realistic we calculate the charging energies E set , E q b and 
E[ nt for the following case: the capacitance of the Josephson junction Cj = 2.0 x 10 _16 F, the capacitances of the normal 
junctions Cn = 1.0 x 10 _17 F and the capacitances of all other capacitors C = 2.5 x 10~ 18 F. We obtain: _E se t ~ 20K, 
E qh « 10K, E iDt w 0.5K. Taking Q qb = 0.35, Q sct = 0.15 and eV = 48K we get AE w 3K, E CB « 14K, and 
V/2 — Ecb ~ 10K. We also assume 2na = 0.1. The measurement time in this regime is t ms w 0.38 x 10 4 ?i/ (ks IK) w 
0.28 x 10~ 7 s. For this choice of parameters we calculate t m ; x numerically, assuming first Ej = 0.1K, and we obtain 
i mix ss 1.4 x 10 5 ?i/(fcB IK) as 1.0 x 10~ 6 s. Thus i m i X /t ms ~ 35 and the separation of peaks should occur much earlier 
than the transitions happen. Indeed, the numerical simulation of the system (p0|-p7|) for those parameters given 
above shows almost ideal separation of peaks (see Fig. ||). Then, we calculate £ m ; x for Ej = 0.25K, and we obtain 
tmix/tms ~ 6. This is a marginal situation. The numerical simulation in this case shows (see Fig. ^) that the peaks, 
first, start to separate, but, later, the valley between the peaks is filled due to the transitions. 

In this paper we have demonstrated that the current through a single-electron transistor can serve as a measurement 
of the quantum state of the q-bit, in the sense that in the case of a superposition of two eigenstates it gives one or the 
other result with the appropriate probabilities. This should be distinguished from another question, namely whether 
it is possible to demonstrate that an eigenstate of a q-bit can actually be a superposition of two different charge 
states, i.e. whether it depends on the mixing angle 77 as described by Eq. (^|). This question has been addressed in the 
experiments of Refs. [ p|,^0|]. They used a setup similar to the one shown in Fig. [j], a single-Cooper-pair box coupled 
to a single-electron transistor. They could demonstrate that the expectation, i.e. the average value of the charge in 
the box varies continuously as a function of the applied gate voltage as follows from (jq) . 

Our theory can also describe the type of measurements performed in Refs. [ |9jjl0|7^ For this purpose we analyze 
the rates in the master equation (19) for general values of the mixing angle 77, relaxing the requirement tan 77 <C 1. 
Then, for our approach to be valid, we must have E- mt -C Ej, so that tane^v <C 1. In this regime each eigenstate of 
the q-bit, |0) or |1), corresponds to a single, though 77-dependent propagation velocity (Tq or Ti). Thus, if the q-bit 
is prepared in one of its eigenstates, then even at the degeneracy point (rj = 7r/2) where the eigenstates are equally 
weighted superpositions of two charge states, one would observe only one peak. We have calculated T as a function 
of 77 using (p8|), (p9|), d30|), and (|||) and obtained curves (not shown here) very similar to those in the experiments. 
It should be added that near degeneracy our setup would not be efficient in projecting onto the eigenstates anymore, 
since the difference between the velocities of the peaks, |To — Til, vanishes near the degeneracy point. 

To conclude we have shown that a single-electron transistor capacitively coupled to a q-bit may serve as a quantum 
measuring device in an accessible range of parameters. We have described the process of measurement by deriving the 
time evolution of the reduced density matrix and we have discussed two dual ways to further reduce it. One way, in 
which the density matrix of the q-bit is obtained, provides the dephasing time, while the other, in which the number 
of tunneled electrons is counted, provides the time of measurement. We have shown that, in our case, the dephasing 
time was shorter than the measurement time, and we have discussed the physical meaning of this result. Finally, we 
have estimated the mixing time, i. e. the time scale on which the transitions induced by the measurement occur. Wc 
have shown that it may be made longer than the measurement time with current technology. 
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FIG. 1. The circuit consisting of a q-bit plus a SET transistor used as a measuring device. 
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FIG. 2. The first order diagram for the transition rates. 
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FIG. 3. P(m,t), the probability that m electrons have tunneled during time t. The parameters are those given in the text, 
Ej — 0.1K. The time is measured in seconds. The initial amplitudes of the q-bit's states: a — V0.75, b = \/0.25. 
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FIG. 4. P(m,t), the probability that m electrons have tunneled during time t. The parameters are those given in the text, 
Ej — 0.25K. The time is measured in seconds. The initial amplitudes of the q-bit's states: a — V0.75, b = V0.25. 
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